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the  atom-target  potential  and  of  force  constant  changes  near  the  surface.  This 
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Introduction 

Previous  theoretical  work  on  one-phonon  processes  has  been  based  on  i  a )  approximate 
force  constants  between  both  bulk  and  surface  target  atoms1'*,  lb)  approximate  atom- 
target  interaction  potentials'* •^•5-6.7  ( usually  a  sum  of  pairwise  potentials9  '*"*  9 )  and  (c)  use 
of  the  distorted  wave  Born  approximation3*' '  •9‘**'9. 

This  work  shows  that  the  creation  (or  annihilation)  of  long-wavelength  phonons  is 
independent  of  (1)  changes  of  force  constants  in  the  surface  region  of  the  unperturbed 
target  and  (2)  the  functional  form  of  the  atom-surface  interaction  potential  V(r\  {u,} ) 
where  r  is  the  position  of  the  colliding  atom  and  iT*  the  displacement  of  the  t‘th  target 
atom.  The  only  requirements  are  that  the  target  force  constants  reproduce  correctly  the 
bulk  elastic  constants  and  that  the  inelastic  scattering  is  aweakv9. 

Therefore  any  disagreement  between  theory  and  experiment  for  long- wavelength  phonons 
cannot  be  attributed  to  inadequacies  of  force  constants  or  the  interaction  potential,  nor 
can  agreement  be  regarded  as  confirmation  of  approximations  (a)  or  (b).  Short-wavelength 
processes  do  provide  tests  of  (a)  and  (b).  Much  of  the  published  theoretical  and  experimen¬ 
tal  work  gives  only  relative  intensities  of  one-phonon  processes  as  functions  of  momentum 
transfer'*’'*''*.  Our  work,  which  gives  exact  results  for  the  intensity  of  long- wavelength  pro¬ 
cesses,  allows  the  determination  of  an  absolute  intensity  scale.  Our  theory  also  shows  that 
the  use,  in  scattering  calculations,  of  an  effective  incident  energy  which  includes  the  well 
depth  is  not  applicable  to  long-wavelength  single-phonon  processes.  This  correction,  intro¬ 
duced  by  Beeby  to  estimate  the  Debye-Waller  factor*9,  has  been  used  also  for  calculating 
the  differential  cross-section*1'*2. 

We  now  present  a  heuristic  explanation  of  why  the  probability  for  the  weak  inelastic 
creation  of  a  long-wavelength  phonon  is  entirely  determined  by  the  bulk  elastic  constants 

of  the  target. 

2 


The  displacements  in  the  surface  region  associated  with  a  (surface  or  bulk)  phonon 
of  small  wave-vector  q  axe  locally  rigid  with  a  coherence  length  of  q~{  both  parallel  and 
perpendicular  to  the  surface.  The  restoring  forces  reside  mainly  in  the  interior  (if  qa  <  1, 
where  a  is  the  lattice  constant)  and  are  determined  by  the  static  bulk  elastic  constants. 
Changes  of  the  force  constants  in  a  few  layers  near  the  surface  axe  irrelevant  to  the  quasi- 
rigid  displacements  in  the  surface  region. 

To  understand  why  the  functional  form  of  V  is  irrelevant  we  consider  (for  simplicity) 
a  surface  with  negligible  corrugation.  The  first-order-distorted-wave-Born  approximation 
matrix  element,  adequate  for  weak  scattering,  is  given  to  lowest  order  in  qa  by 

y.  .  .  *(t.  tnii*.  )  (1) 

iUK+Q.E±h~-.K.E  y'K.E'  dz  K,E) 

provided  that  qa  <  1.  Here  V'(c)  =  V'(r:  {0}),  Q  and  u  are  the  parallel  momentum  and 
frequency  of  the  phonon,  and  K  and  E  are  the  parallel  momentum  and  energy  of  the 
incident  particle.  If  $  g  £  is  normalized  to  unit  incident  perpendicular  flux,  the  right  hand 
side  represents  the  total  force  exerted  by  the  rigid  target  per  unit  time  and  area.  This  is 
rigorously  given  by  the  change  of  perpendicular  momentum  per  unit  time  and  area,  i.e. 

(**.*•  (2) 
independent  of  the  functional  form  of  V(z). 

Our  results  are  exact  in  the  limit  where  the  following  conditions  are  satisfied  (in  addi¬ 
tion  to  qa  <£  1): 

JL«!  (3a)  BfJL)  «1  (36)  £<»~'D  (3c)  kBT<huD  (3d) 

M~b  M  \huDJ 

Here  c  is  a  characteristic  velocity  of  a  phonon  in  the  target,  b  is  the  range  ot  the  atom- 
target  potential13  and  V0  is  its  depth,  m  is  the  mass  of  the  colliding  atom,  M  is  the  mass 
of  a  target  atom,  is  the  Debye  frequency,  and  T  is  the  temperature  of  the  target.1 
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The  condition  (3a)  allows  one  to  stop  the  expansion  of  V  after  the  first  order  term 

V(r;  {ut})  =  V(z)  +  Uj  ■  {0})  , 


(4) 


Condition  (3b)  ensures  that  in  calculating  the  inelastic  scattering,  the  last  term  in  (4)  can 
be  treated  to  lowest  order  (Fermi’s  golden  rule).  The  equations  (3a)  and  (3b)  are  derived 
in  detail  in  the  appendix. 


Derivation  of  the  Probability  for  Production  of  Long- Wavelength  Phonons 


To  simplify  the  presentation,  we  will  make  the  following  assumptions  in  this  section 
(see  concluding  remarks):  (1)  negligible  diffraction  scattering;  (2)  isotropy  of  the  target; 
(3)  interaction  only  with  the  uppermost  layer  of  target  atoms.  The  units  are  chosen  so 
that  the  volume  of  the  target  and  the  area  of  the  interacting  surface  of  the  target  are  1. 

To  calculate  the  one-phonon  amplitude  we  divide  the  Hamiltonian  into  an  unperturbed 
Hamiltonian  H0,  which  ignores  phonon  coupling  to  the  colliding  atom,  and  which 
couples  the  atom  linearly  to  the  phonons**. 


H  H particle  ~b  H target  "b  L  (r;  {  u,  }  )  —  Hq  +  H i  -f-  0(  U2)  (5) 


where 


2  * 

=  2m  +  +  W-h  #1  =  X]  gU(  ~^Ui<a 


i, a 


(6) 


Here  a  denotes  a  cartesian  coordinate,  q  is  the  wavevector  of  a  phonon  and  v  is  its  branch. 


The  colliding  atom  wavefunctions  in  the  potential  V(c)  have  the  form  ^ £  ^(r)  = 

r**  p  • 

em  n<^K  £■(*)  ^or  atom  energy  E.  K  is  momentum  parallel  to  the  surface,  and  R  is 
displacement  parallel  to  the  surface.  ct>^  E(z)  satisfies  the  equation 


h2d2 

2m 


+  V(z)-(E- 


h2I<2 

2m 


) 


^K,E^  ~  °' 


(7) 
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tig  £{z)  has  the  asymptotic  form  (m/hk;) 2  e~,k'z  +  A0e,k‘z  as  c  -»  3C.  where  Jt?  +  K2  = 
2  mE/h2. 

Before  calculation,  we  expand  the  ut-  in  normal  modes, 

&\r 

=  E  +h  c.  (8) 

i.a.q,i/  1,0 

which,  for  convenience,  we  rewrite  as 


g,£/.a  i  ,LZ 


-tq(R,-R) 


+h.c. 


If  the  wavelength  of  the  phonon  is  much  larger  than  the  range  of  the  coupling  potential 
(similar  to  a  and  b),  dV /du(  vanishes  unless  q(Ri  —  R)  <£  1  and  the  second  exponential 
can  be  ignored.  This  yields 

H'  =  e~iq '%>.<*  H  +h-c-  (10) 


Since  V(r,  {£,})  depends  only  on  r  -  iT, , 

o-y^Tj—R' 

drQ  ~  5ui,a 

Since  ( :)/di  =  dV(z)/dy  =  0,  simplifies  to 

rr  -tafidV,  , 


=  +h-c. 


The  differential  probability  of  scattering  into  final  angle  Q  with  energy  E\  per  unit 
incident  perpendicular  flux  is  given  by  the  golden  rule: 


d2P  2; r 
dQdE'  ~  h 


=  f  | (. . .  n,-„  ±  1 . . .  |  Jv  VR,  £,Hx* R  Edr\ . . .  n,y . .  -)l 2p(E')  x 


where  is  the  Bose  occupation  factor  for  the  phonon  q,  v,  and  p(E')  is  the  density  per 
unit  solid  angle  of  states  ’P g,  ^ , .  p{E')  is  computed  with  continuum  normalization  in  the 
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r-direction  but  box  normalization  in  the  x  —  y  plane.  The  matrix  element  for  E'  <  E 


(creation),  appearing  in  equation  (13).  is 


i  -  „  f  dV 

nqu  ~  q  fct'q.i'.z  J  £i  ~q~®K  E^Z‘ 


where  u$M  =  +  a_q- „) 

The  difference  in  energy  and  parallel  momentum  between  the  initial  and  final  colliding 
atom  states  is  very  small.  As  a  result,  we  may  replace  0*r.t  by  0*-  .  The  integral  is 

now  ( dV/dz ),  the  expectation  value  of  the  force.  It  is  elementary  to  show,  and  it  has  been 
remarked  in  the  context  of  other  surface  scattering  problems,15  that 

=  M  +  1  \d:oml:)ds<t>(:)-(&0m(zmz)]  = -2  hkz.  (15) 

'  dz  at  2m  L  r=oo 

Using  (14)  and  (15)  in  (13)  results  in 

-  £';<?)<",->  + 1}.  lie) 

g(E  —  E’\Q)  is  the  density  per  unit  energy  of  phonons  with  energy  hu  =  E  -  E1  and 
surface  wavevector  Q. 

The  value  of  U~u  .  in  the  long- wavelength  limit  for  an  isotropic  target,  like  tungsten, 

is  well-known7.  There  are  three  types  of  phonons  with  displacement  in  the  z-direction: 

(1)  Rayleigh  waves16,  with  a  dispersion  relation  =  £ ctQ  (c<  is  the  transverse  speed  of 

sound.  Q  is  the  momentum  of  the  Rayleigh  wave,  and  0  <  £  <  1);  (2)  mixed  phonons  ■  , 

with  dispersion  relations  ctQ  <  ^  <  ceQ  {ct  is  the  longitudinal  speed  of  sound);  (3)  bulk 

waves1'  18  with  dispersion  relations  ^  >  c^Q. 

For  Rayleigh  waves^,  Uq  —  hF(a)  /(pp)i,  where  p  is  the  density  of  the  target,  F 

is  a  dimensionless  number  depending  only  on  the  Poisson  ratio  <r,  and  \,p  are  the  Lame 

on 

coefficients.  For  tungsten  F  has  the  value  0.160.  For  a  general  isotropic  target  , 


2  -  £2 


2-^ 


(2-£2)2(2-(£q)2) 

8(l-(£a)2)f 
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where 


/T  — 2  a 

2-2<7 


and  in  the  interval  0  <  f  <  1.  satisfies 


^  -  S£4  +  S£2(3  -  2a2)  -  16(1  -  a2)  =  0.  (19) 

There  is  only  one  Rayleigh  wave  for  each  value  of  the  surface  momentum,  so  g(hus\  Q)  = 
6(h(uj  -  (ctQ)).  When  the  integration  over  E'  is  performed,  we  find 

dP  hkhF  , 

5a  =  ',,'  .("fo  +  D-  (20) 

au  -‘(mp)2 

For  mixed-mode  phonons.  C  =  hG  (&,  3)  /  pu,'  where  3  =  -jjfc^Q.  G  is  a  dimensionless 
number  that  depends  on  the  Poisson  ratio  and  the  ratio,  3.  of  the  “surface  velocity"  jj/Q 
toe*. 

r  4J'2f  J2  -  1K1  -(3a?) 

[16(J2  -  1)(1  _(Jq)2)  +  (2  -  J2)4]  •  t21) 

Since  only  the  transverse  part  of  a  mixed-mode  phonon  extends  through  the  bulk,  the 
mixed  mode  density  of  states  is  identical  to  that  of  the  transverse  modes, 

g(hu;:Q)  = - ^ - r-  (22) 

xhMJi  -  1)}* 

For  bulk  phonons,  the  only  difference  in  U~,  from  mixed  modes  is  that 

r  t(jq)2-l}J2  1 

G~  U(a2-l)  +  jwJ'  (-3) 

Finally,  the  density  of  states  for  bulk  modes  is 

A\  y/P  i  d5a2  -  4J(1  -  a2)  1 

g(tv\Q)  =  <  -==  —  r  ;  ==~  =Z== - — (24) 


V?  \  -  1  [WP  -W(0<*)2  -  1  +  (2  -  32)2 


Experimental  Implications 


We  present  below  three  graphs  (Figures  1-3)  of  dP/dQdE'  for  17  meV  Helium  atoms 
scattering  off  of  a  70A'  tungsten  surface21  ( kgT  ~  hug/ 4)  which  happens  to  be  weakly 
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corrugated.  We  chose  tungsten  because  it  is  very  isotropic22  and  the  mass  ratio  m/M  <  1. 
The  results  for  Rayleigh  phonons  are  presented  on  the  same  graphs  by  displaying  dP/dCl, 

Concluding  Remarks 

Our  results  were  derived  under  three  simplifying  assumptions.  The  results  remain 
unchanged  if  the  third  assumption  (interaction  with  only  the  top  layer)  is  dropped.  For 
systems  to  which  the  first  two  assumptions  (negligible  diffraction,  target  isotropy)  do  not 
apply,  exact  model- independent  results  of  the  same  general  nature  are  still  obtained  but 
have  a  more  complex  form.  The  more  general  results  will  be  presented  in  a  separate 
publication23. 

Appendix 

An  estimate  of  higher-order  contributions  to  the  one-phonon  amplitude  provides  re¬ 
assurance  that  the  first-order  Born  contribution  dominates.  The  two  types  of  new  terms 
come  from  1)  higher  order  perturbation  of  the  linear  coupling  term  H\  and  2)  expansion 
of  the  potential  V(r;  {if})  to  higher  order  in  the  phonon  coordinates  {tf}.  For  both  of 
these  the  next  contributing  term  is  third  order.  We  find  two  parameters,  A  and  a.  which 
characterise  the  higher  order  terms  of  type  (1)  and  (2)  respectively.  These  parameters 
are  estimates  of  the  ratio  of  the  third-order  processes  to  the  first-order  processes.  The 
requirement  a  <C  1  is  equation  (3a)  and  A  <  1  is  equation  (3b). 

Consider  the  third-order  processes  of  type  (1).  These  processes  require  a  transition  to 
an  intermediate  state  with  a  second  phonon  and  then  the  destruction  of  the  same  phonon. 
(Processes  which  interchange  the  creation  and  annihilation  events  are  characterized  by 
the  same  parameter  A).  The  transition  to  an  intermediate  state  contributes  a  factor  of 
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(H i)S/AE',  here  (H i)  is  a  characteristic  matrix  element  of  the  linear  coupling  term  in 
the  Hamiltonian,  A E  is  a  characteristic  energy  denominator,  and  .V  is  the  number  of 
intermediate  states  in  the  characteristic  range  of  energies.  If  the  intermediate  particle 
state  is  unbound  then  .V  ~  p(otai{E)SE\  here  Pt0tal{E)  is  the  density  of  states  of  the 
atom-target  system  at  a  characteristic  energy.  Since  the  destruction  event  allows  no  choice 
of  the  target's  final  state,  the  characteristic  ratio  of  the  third  order  processes  to  the  first- 
order  processes,  A.  is  (Hx)' ptotapE)ppariiclli(EpartXcu)-  The  overline  in  equations  and  the 
adjective  'characteristic-  in  the  text  will  henceforth  be  dropped. 


J  Pparticle^  ^p)Pphonon(^  ~  Ep)dEp  —  fic‘2yy~£^  huijp 


The  density  of  states  Ptotal(E)  is 
•E 
1 0 

fix)  is  bounded  from  above  by  a0x~  for  x  1  and  by  a0 ©  otherwise;  aQ  and  a0 o  are 
constants  of  order  unity.  The  maximum  value  of  Pt0tal(E)  thus  occurs  for  E  ~  ftw£>.  We 
will  use  Ptotal(h~D)  ‘n  our  calculation  of  A. 


The  matrix  element  (H i)  is 

[  ,  ,dV(r) 

/  dri/q  — — — l  d'2  ~ 


VQ 


HAC 

M^D 


(A2) 


U  is  a  phonon  amplitude  and  .4c  is  the  area  of  the  surface  unit  cell.  Characteristic  phonons 
will  couple  particle  states  that  differ  by  a  momentum  of  1/a.  Intermediate  particle  energies 
will  thus  be  of  order  £2/2ma2  ~  h^'D- 

Combining  these  factors  yields 


A  ~  (Hi)"2 Ptotali^ D)Pparticle(^lu! D) 


m  Vq 

a/(^d)2‘ 


(-43) 


_ o 

For  helium  on  tungsten  A  «  7  •  10  . 

The  contribution  from  bulk  phonons  to  A  is  comparable  to  the  contribution  from  surface 
phonons.  For  each  surface  momentum  Q  there  is  one  surface  phonon  and  there  are  Nt 
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bulk  phonons.  Nz  is  the  number  of  layers  of  the  solid  in  the  z  direction.  However,  the 

WO 

additional  factor  l/.\y  in  Uq  cancels  this  Sz  dependence. 

The  third-order  processes  whose  intermediate  particle  states  are  bound  are  also  charac¬ 
terized  by  A.  For  this  process  (H\)  ~  Y0Ua~1/2  and  ptotiE)  ~  m/h2  so  A  ~ 

AI  ( riuj  hiuj  q 

Here  £3  is  the  energy  of  the  bound  state.  There  are  no  bound  stafes  with  £3  »  hu;D, 
Other  third-order  processes  not  explicitly  discussed  here  are  quickly  seen,  by  similar  argu¬ 
mentation.  to  be  characterized  by  A. 

The  other  parameter  in  the  expansion  is  the  ratio  of  the  u3  terms  to  the  u  terms 
in  the  expansion  of  V(r ;  {u}>.  both  types  of  terms  calculated  in  first-order  perturbation 
theory.  This  ratio  is  (cfiV/d hi)U2 / (dV/ dz) .  If  the  potential  V(r)  has  range  b  then  a  ~ 
(V/l)2.  For  a  surface  phonon  of  the  Debye  frequency,  a  ~  ( h/2Mcb ).  The  previous 
arguments  concerning  relative  contributions  of  bulk  and  surface  phonons  apply.  For  helium 
on  tungsten,  a  %  10-4. 
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Figure  Captions 


1.  Graph  of  d?P/dCldE'  (dP/d?l  for  interaction  with  Rayleigh  phonons)  for  17  meV  Helium 
atoms  scattered  by  a  70A'  tungsten  surface.  The  incident  angle  of  the  beam  is  44°  from 
the  perpendicular.  The  final  angle  is  46.5°  and  is  in  the  sagittal  plane.  The  arrows  point 
to  the  peaks  which  axe  caused  by  Rayleigh  phonons. 

2.  The  same  parameters  are  used  as  in  Fig.  1.  except  that  the  final  angle  is  47°. 

3.  The  same  parameters  are  used  as  Fig.  2.  except  that  the  outgoing  beam  is  displaced 
from  the  sagittal  plane  by  0.65°. 
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